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1. Weak radiative B-meson decay is known to be a very sensitive probe of new physics |]T], Q]. 
Heavy Quark Effective Theory tells us that inclusive 5-meson decay rate into charmless 
hadrons and the photon is well approximated by the corresponding partonic decay rate 

r(5^X,7) ~r(6^X,7). (1) 

The accuracy of this approximation is expected to be better than 10% 0. 

The inclusive branching ratio B{b^Xs'~f) was extracted two years ago from CLEO mea- 
surements of weak radiative 5-meson decay. It amounts to 

B{b^Xs-f) = (2.32 ± 0.57 ± 0.35) x 10"^ (2) 

In the forthcoming few years, much more precise measurements of B{B^Xs'~f) are expected 
from the upgraded CLEO detector, as well as from the B-factories at SLAG and KEK. Ac- 
quiring experimental accuracy of below 10% is conceivable. Thus, the goal on the theoretical 
side is to reach the same accuracy in perturbative calculations of b^Xs'y decay rate. This 
requires performing a complete calculation at the next-to-leading order (NLO) in QCD |l], |^. 
All the next-to-leading contributions to b^Xs'y are collected for the first time in the present 
letter .0 

One of the most complex ingredients of the NLO calculation consists in finding three- 
loop anomalous dimensions in the effective theory used for resummation of large logarithms 
\n{M^/ml). This is the main new result presented in this letter. However, we have recal- 
culated all the previously found anomalous dimensions, too. 

The remaining ingredients of the NLO calculation are already present in the literature 
0^113 • published results for the two-loop matrix elements calculated by Greub, 



Hurth and Wyler Bremsstrahlung corrections obtained by Ali, Greub [jT5|, |T6[ and Pott 



17| , and the matching conditions found by Adel and Yao ||13|| . 

The following three sections of the present paper are devoted to presenting the complete 
NLO formulae for b^Xg'y. Next, we analyze them numerically. In the end, we discuss 
predictions for the branching ratio B{B^Xs'y) including a discussion of the relevant uncer- 
tainties. 



^ except for the unknown but negligible two-loop matrix elements of the penguin operators which we 
denote further by P3,...,Pq. 
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2. The analysis of b^Xs'j decay begins with introducing an effective hamiltonian 

AC ^ 

Heff = ~^v*v,,J:c^i^^)P^i^^) (3) 

where Vij are elements of the CKM matrix, are the relevant operators and Cj(yu) are 

the corresponding Wilson coefficients. The complete set of physical operators necessary for 
6— >Xs7 decay is the following: 



Pi = 


{sl1^,T^cl){cl1^T%l) 


P2 = 


{sLl^iCL){cLl^hL) 


P3 = 




P,= 




P,= 




P6 = 




P7 = 




P8 = 





(4) 



where T" stand for SU (3)coior generators. The small CKM matrix element Vub as well as the 
s-quark mass are neglected in the present paper. 

Our basis of four-quark operators is somewhat different than the standard one used in 
refs. Jl], H, |l8l, although the two bases are physically equivalent. We would like to stress 
that in our basis, Dirac traces containing 75 do not arise in effective theory calculations 
performed at the leading order in the Fermi coupling Gp (but to all orders in QCD). This 
allows us to consistently use fully anticommuting 75 in dimensional regularization, which 
greatly simplifies multiloop calculations. Our choice of color structures is dictated only by 
convenience in computer algebra applications. 

Throughout the paper, we use the MS scheme with fully anticommuting 75. Such a 
scheme is not uniquely defined until one chooses a specific form for the so-called evanescent 
operators which algebraically vanish in four dimensions ^ . In our three- loop calculation, 
eight evanescent operators were necessary. We list them explicitly in Appendix A. 

Resummation of large logarithms ln(M^/m^) is achieved by evolving the coefficients 
from /i = Mw to fi = — f^b according to the Renormalization Group Equations 
(RGE). Instead of the original coefficients Ci{fi), it is convenient to use certain linear com- 
binations of them, called the "effective coefficients" ^ Q 

{Ci{fi) forz = l,...,6 

C7i^i) + EliyiQ{ix) iort = 7 (5) 
Cs{^i) + Etl^^C,{^I) for 2 = 8. 
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The numbers yi and Zi are defined so that the leading-order h ^ s'~) and h ^ s gluon 
matrix elements of the effective hamiltonian are proportional to the leading-order terms in 
Cj^^ and Cg-^-^, respectively In the MS scheme with fully anticommuting 75, we have 
y = (0, 0, -i, -|, -f , -f ) and z = (0, 0, 1, -i, 20, The leading-order contributions 

to the effective coefficients Cl^^{fi) are regularization- and renormalization-scheme indepen- 
dent, which would not be true for the original coefficients Ct{^) and Csifi). 
The effective coefficients evolve according to their RGE 



(6) 



driven by the anomalous dimension matrix ^^^^ (/i). One expands this matrix perturbatively 
as follows: 



(7) 



An ' ' (47r)2 

The matrix ^^^^^ff is renormalization-scheme independent, while ^^^'^^^f is not. In the MS 
scheme with fully anticommuting 75 (and with our choice of evanescent operators) we obtain 



and 



7 









-4 


8 
3 





2 
9 





208 
243 


173 " 
162 










12 








4 
3 





416 
81 


70 
27 



















52 
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2 


176 
81 


14 
27 














- 



40 


100 


4 5 


152 


587 




^(0)e// 






9 




9 
256 
3 


9 6 

20 


243 
6272 
81 


162 

6596 
27 
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256 


56 


40 2 


4624 


4772 










9 


9 


9 3 


243 


81 

























32 
3 




























32 
9 


28 
3 




355 


502 




1412 




1369 


134 


35 


818 


3779 


9 


27 




243 




243 


243 


162 


243 


324 


35 


28 




416 




1280 


56 


35 


508 


1841 


3 




3 




81 




81 


81 


27 


81 


108 










4468 




31469 


400 


3373 


22348 


10178 




81 




81 


81 


108 


243 


81 










8158 
243 




59399 
243 


269 
486 


12899 
648 


17584 
243 


172471 
648 










251680 




128648 


23836 


6106 


1183696 


2901296 




81 




81 


81 


27 


729 


243 










58640 




26348 


14324 


2551 


2480344 


3296257 




243 




243 


243 


162 


2187 


729 
























4688 
27 



























2192 
81 


4063 
27 



(8) 



(9) 



^ These numbers are different than in section 4 of ref. [Q because we use a different basis of four-quark 
operators here. 
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The matrix the main new result of this paper. Its evaluation required performing 

three-loop renormalization of the effective theory (|^). Details of this calculation will be given 
in a forthcoming publication . 



We expand the coefficients in powers of as follows: 



(10) 



At = My/i the values of the coefficients are found by matching the effective theory ampli- 
tudes with the full Standard Model ones. The well-known leading-order results read |jl|, ^, ^ 





1 

■ix^-2x^ 1 _l_ -8a;^-5a;^+7x 
4(x-l)4 ^^^-^ 24(x-l)3 

1)4 "^-^"T 8(a;-l)3 



k 4(1-1) 



for i = 1, 3, 4, 5, 6 
for 1 = 2 
for i = 7 

for z = 8, 



(11) 



where 



X 



^t,pole 



1 



3 TT 



, , , >^ o (12) 

The next-to-leading contributions to the four-quark operator coefficients at /i = M\y can 
be found either by transforming the results of ref. to our operator basis, or by a direct 
computation. We find 

for i = 1 

i = : 

for i = 4, 
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E{x)- 



for 2 = 2, 3, 5, 6 



(13) 



where 



, x l8-llx-x2 x2 15- 16x + 4x2 , 2, 

E(x) = — -T — - H ^ , ^Inx Inx. 

^ ' 12(1 -a;)'^ 6(1 -x)4 3 



(14) 



The next-to-leading contributions to Cj^^ {Mw) and C^^^ {M\y) can be extracted from 



ref. according to the following equations: 



TT 



C^^'f^Mw) = + Fix) 



As--{E(x)-- 



(15) 



(16) 



where the quantities Cq^^ and Cq^^ are given in eqn. (5) of ref. [|T^ with fi = M^r. The terms 
denoted by A7 and Ag stand for shifts one needs to make when passing from the so-called 



R* renormalization scheme used in ref. to the MS scheme used here 



A, 



47r 

a. 



Cf\Mw){x ^x- —x\nx) - Cr{Mw){x) 

71 



2a. 



for i = 7, 



(17) 



Only the leading (zeroth-order) term in as needs to be retained in the above equation. 
The resulting explicit expressions for Cj^^{Mw) and C^^^{Mw) in the MS scheme read 
,(i)ej/ _ -16x^ - 122x3 + 80x2 - 8x . 6x^ + 46x^-28x2 ^ 



= 9(^' 3(x-l)^ '"-^ 

-102x5 - 588x^ - 2262x3 + 3244x2 - 1364x + 208 
81(x — 1)^ 

1646x^ + 12205x3 - 10740x2 + 2509x - 436 



486(x - 1)4 

_ -4x4 + 40x3 + 41x2 + X . -17x3-31x2 2 



" f^''") - ['--.)^ 2(.-ip ^ 

-210x5 + 1086x4 + 4893x3 + 2857x2 - 1994x + 280 , 

H ; Inx 

216(x- 1)5 

737x4 _ 14102x3 - 28209x2 + 610x - 508 

^ 1296(x- 1)4 ^^^^ 

Having specified the initial conditions at /i = M^y, we are ready to write the solution to 
the Renormahzation Group Equations (|) for C,^^^(/i,) 

= V^C?\Mw) + I - r/i) ci'\Mw) + E (20) 



/297664 16 7164416 m 256868 37 6698884 39\^(o).,. x 

+ 7723 77 23 H 77 23 r723 Cs (Mm/) 

V 14283 ' 357075 ' 14283 ' 357075 ' / ^ ^ ^ 



37208 - - ^ 



4761 



(ryi - r^i) (Mh/) + + + ^?*^)^"% (21) 



i=l 

where rj = as{M\y) /agifJ^b), and the numbers Oj-Zij follows 



Cti = ( 


14 
23' 


16 
23' 


6 

23' 


12 
23' 


0.4086, 


-0.4230, 


-0.8994, 


0.1456) 


ei = ( 


4661194 
816831 ' 


8516 
2217' 


0, 


0, 


-1.9043, 


-0.1008, 


0.1216, 


0.0183) 


/i = ( 


-17.3023, 


8.5027, 


4.5508, 


0.7519, 


2.0040, 


0.7476, 


-0.5385, 


0.0914) 


= ( 


14.8088, 


-10.8090, 


-0.8740, 


0.4218, 


-2.9347, 


0.3971, 


0.1600, 


0.0225) 


hi = ( 


626126 
272277' 


56281 
51730 ' 


3 

7' 


1 

14' 


-0.6494, 


-0.0380, 


-0.0186, 


-0.0057) 



As far as the other Wilson coefficients are concerned, only the leading-order contributions 
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^0.1456 
^0.1456 
„0.1456 



to them are necessary in the complete NLO analysis of b^Xg'y. We obtain^ 

I \ 1 -12 9 A 

M = 3^ 23 + -7^23, 

{^^^) = -^r^-i + ^r^i-O.OeSV'^os^ + 0.0595r7-°-^230 _ o.02 1 87y-°-^9^^ + O.OSSSr^^ 
(/ife) = 1^"^ + +0.0237r/0-4°s6 _ o.0173r/-o-4230 _ o.l336r/-o-8994 _ o.0316r/° 

(^^) = Y^T^-i - ^T]^ +0.009V-^os*5 - 0.0100r7-°-4230 ^ o.OOlOr/-^-^^^^ - 0.0017r/° 
(^j,) = -^r/-i - ^r/A +0.0108r7°-^°^6 + 0.0163r7-°-4230 ^ o.0103r/-°-^99^ + 0.0023r70-i^5^ 

cr-(..,^(cr(M.). Ill) 

- 0.9135r7™ + 0.0873r/-°-^230 _ o.0571r7-o-899^ + 0.0209r/°-i^5^(22) 

In our numerical analysis, the values of as{fi) in all the above formulae are calculated 
with use of the NLO expression for the strong coupling constant 



An v{iJ,) 



(23) 



where 



= 1 -/5o^^ In (^—j, (24) 

/3o = f and pi = ^. 

As easily can be seen from eqns. (p!0D-(p4D, the Wilson coefficients at ^ = fit, depend on 
only five parameters: Mz, Mw, Ois{Mz), rrit^poie and /if,. In our numerical analysis, the first 



two of them are fixed to Mz = 91.187 GeV and Mw = 80.33 GeV [^. For the remaining 
three, we take 

as{Mz) = 0.118 ±0.003 [24], (25) 
mt,poie = 175 ±6 GeV ^ (26) 

fib e [2.5 GeV, 10 GeV]. (27) 
For /if, = 5 GeV and central values of the other parameters, we find as{fib) = 0.212, 

C^°^'^^(/ib) = (-0.480, 1.023, -0.0045, -0.0640, 0.0004, 0.0009, -0.312, -0.148)^ 
Ci^^^^il^b) = -0.995 + 1.443 = 0.448. (2^ 



In the intermediate step of the latter equation, we have shown relative importance of two con- 

(w) 



tributions to C^^^'^^^Hb)- The first of them originates from terms proportional to Cj^^^^\M]a 



Analogous expressions in ref. Q are somewhat different, because we use a different basis of four-quark 
operators fiere. 
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in eqn. (|2T|), i.e. it is due to the NLO matching. The second of them comes 
from the remaining terms, i.e. it is driven by the next-to-leading anomalous dimensions. 
Quite accidentally, the two terms tend to cancel each other. In effect, the total influence 
of Cj^^^^\fib) on the b—'-Xsj decay rate does not exceed 6%. (Explicit expressions for the 
decay rate are given below). However, if the relative sign between the two contributions was 
positive, then the decay rate would be affected by almost 30%. This observation illustrates 
the importance of explicit calculation of both the NLO matching conditions and the NLO 
anomalous dimensions. 

3. At this point, we are ready to use the calculated coefficients at yU = /i^ ~ mf, as an input in 
the NLO expressions for b—>-Xs'~f decay rate given by Greub, Hurth and Wyler 0. However, 
we make an explicit lower cut on the photon energy in the Bremsstrahlung correction 

e^>(i-5)e::;""-' = (i-5)^. (29) 

We write the decay rate as follows: 

r[6^X,7]^->(i-'^)^r^ = r[6^s7] + r[b s-f gluonf->^'-'^''r'' = 

\V*V,b\'ml^,MMsi^>>) (\DnA), (30) 



327r4 



where 



An 

and 



i=l 



, (0)e// , "^6 

+ Ti7 In — 



(31) 



i<j 

The contribution denoted by is independent of the cutoff parameter 6. By differenti- 
ating D with respect to /if, and using the RGE (^, one can easily find out that the leading 
/ifc-dependence cancels out in D. When calculating |Z)p in our numerical analysis, we con- 
sistently set the O^a^) term to zero. 

The contribution denoted by A vanishes in the limit a^-^O. The first term in A contains 
the exponentiated (infrared) logarithms of 6 which remain after the IR divergences cancel 
between the virtual and Bremsstrahlung corrections to b-^Xg'y. 

The terms denoted by rj in D depend on what convention is chosen for additive constants 
in the functions fijiS) and As (5)- We fix this convention by requiring that all fij{S) vanish 
in the formal limit 6^0. 
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In order to complete the full presentation of b^Xg'y results, it remains to give explicitly 
the constants and the functions fij{S). 

The constants rg and r2 are exactly as given in ref. 



r2 



4 

- — (-33 



27r^ - Gin), 



(33) 



(-833 + IUtt^z^/^ 

243 ^ 



+ 
+ 



1728 - 180n^ - 1296C(3) + (1296 - 3247r2)L + 108L^ + 36L" 



IGvri f 



648 + 727r2 + (432 - 2167r2)L + 36L^ 

z 



-54 - SAtt^ + 1092L - 756L2 



45 - 37T^ + 9L + 9L2 



+ 



-37r^ + 9L' 



z 



+ [28-12L]z^} + 0{z^L^), 



where z = ml /ml and L = lnz. In our operator basis 
but is proportional to r2 



ri 



1 
6 



r2- 



The constant ri turns out to be^ 



10 8 2 

rv = IT . 

3 9 



(34) 

the constant ri does not vanish, 

(35) 

(36) 



The quantities r^, rg originate from two-loop matrix elements of the penguin operators 
P3, Pg. They remain the only still unknown elements of the formally complete NLO 
analysis of b—>-Xs'~f. However, the Wilson coefficients Ci are very small for i = 3, 6. In 
consequence, setting rs, rg to zero (which we do in our numerical analysis) has a negligible 
effect, i.e. it can potentially affect the decay rate by only about 1%. 

Let us now turn to the functions fij{6). One can easily obtain frriS) from eqns. (23) 
and (24) of ref. [0. (These equations were also the source for our rj and the exponent in 
eqn. (B^) above.) 



/77(5) = y5 + ^5^-^5^ 



-6{6-4) In 6. 



(37) 



The function fssi^) can be found by integrating eqn. (18) of the same paper 
fssiS) = i-(-21n— [52 + 25 + 41n(l-5) 



27 



^ It differs from the one given in ref. 
to /77(1) was absorbed into r^. 



where a formal hmit 5 — > 1 was taken and the contribution due 



+4Li2(l - 5) - ^ - 5(2 + 5) ln5 + 8 ln(l -5)- 



35' + 76 



{31 



The logarithm of in the above equation is the only point in our analysis, at which the 
s-quark mass cannot be neglected. Collinear divergences which arise in the massless s-quark 
limit have been discussed in ref. The properly resumed photon spectrum was there found 
to be finite in the limit, and suppressed with respect to the case when \n{mh/ms) 

in fssi^) is used. Even with as small as 0.1 GeV, the term proportional to \n{mi,/ms) 
is negligible, i.e. it affects the decay rate by less than 1%. Therefore, for simplicity, we will 
just use the naive expression (^) with mf,/ms = 50 in our numerical analysis. 

The functions /22(5), /27(5), /28(5) and fjsiS) can be found from appendix B of ref. 
After performing some of the phase-space integrations, one finds 



/22(5) 



IQz 



5 1 dt{l- zt) 





Git) 1 
t 2 



i-l/z 

+ dt(l- ztf 

J{l~S)/z 
1/z 



Gjt) 1 
t 2 

t 



6 1 dtRei G(t) + -] + r dt(l- zt)Re ( G(t) + - 
V 2/ J(i-s)/z \ 2 



-3 /27(5), 
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TT^ . 9. 1.2 , 1 ,3 



Li2(l -6)---6\n6 + -6--6' + -6 



(39) 



where, as before, z = ml /ml and 
G{t) = 



-2 arctan' ^t/{A-t), for t < 4 

-7rV2 + 2 ln'[(v^ + Vt^)/2] - 2m ln[(v^ + Vt^)/2], for t > 4. 



(40) 



In our operator basis @, the functions fij{S) do not vanish, but are proportional to the 
functions /2j(5) 

/ii(5) = ^22(6), fuiS) = -1/22(5), /i7(5) = -1/27(5), /l8(5) = -|/28(5). 



The only functions fij{S) which have not yet been given explicitly are the ones with at 



least one index corresponding to the penguin operators P3, Pg- As checked in ref. ||17|| , 
these functions cannot affect the decay rate by more than 2%. Therefore, we neglect them 
in our numerical analysis. However, for completeness, a formula from which they can be 
obtained is given in our Appendix B. 



4. The decay rate given in eqn. (|30| ) suffers from large uncertainties due to 



m 



b,pole 



and the 



CKM angles. One can get rid of them by normalizing b^X^'y decay rate to the semileptonic 



decay rate of the b-quark 



where 



(7(2) = l-8z + 8^=^-/-122Mn;2 ["^ = I^Ie^^"] (42) 



is the phase-space factor, and 



= 1 - (43) 
3vr 5f(z) 



is a sizable next-to- leading QCD correction to the semileptonic decay The function 

h{z) has been given analytically in ref. ^8 



/25 239 25 ,\ , / 4 , 17 o\ o, 2 . 9 

/i(;2) = -(1 - z^) i — - —z + —z^j + ;z In ;z f 20 + 90z - -z^ + —z^j + z^ In^ z (36 + z^ 

+ (1 - z^) (—-—z + —zA ln(l -z)~ 4(1 + 30^2 + z"") In z ln(l - z) 



3 3 3 

:i + 16^2 + Z^)[6U2{Z) - TT^] - 32^3/2(1 + z) 



7r2-4Li2(v^)+4Li2(-v^)-21nz In f^^^ 



Thus, the final perturbative quantity we consider is the ratio 

itq«arfc((>) = ^^77 7^— 1 = i.. |2 TT^ ( + A) • (44) 

where D and A are given in eqns. (|3lD and (0), respectively, and 



1 mfc(/x = mb) y _ 1 / _ 8 ag(mfc) ' 



(45) 



5. Let us now turn to the numerical results. Besides the five parameters listed above 
eqn. (|25|) , the quantity RquarkiS) in eqn. (^^ depends on a few more SM parameters. They 
are the following: 

(i) The ratio mc^poiel'mh,-poie- Similarly to ref. 0, we will use 

0.29 ±0.02, (46) 



^b,pole 

which is obtained from m},^poie = 4.8 ± 0.15 GeV and rrih^poie — t^c,poie = 3.40 GeV. 

(ii) The 6-quark mass itself. Apart from the ratio mc/mb, the 6-quark mass enters our 
formulae only in two places: in the explicit logarithm in D ( [31]) and in the factor F (^5]), as 
the argument of ag- In both cases, the mt-dependent terms are next-to- leading, and their 
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mfe-dependence is logarithmic. Changing rrih from 4.6 GeV to 5 GeV in these places affects 
b^Xs'y decay rate by less than 2%. Thus, it does not really matter whether one uses the 
pole mass or the MS mass there. We choose to use rrib^poie = 4.8 ± 0.15 GeV. 

(iii) The electromagnetic coupling constant Oem- It is not a priori known whether this 
constant should be renormalized at /i ~ or /i ~ Mw, because the QED corrections have 
not been included. We allow aem to vary between aemijnh) and aem(^iy), i-e. we take 



130.3 ±2.3. 



(iv) The ratio of CKM factors. We will use 



Vr 



cb 



0.976 ±0.010 



(47) 



given in ref. [|29|. It corresponds to gaussian error analysis, which is in accordance with 
interpreting the quoted error as a la uncertainty. 

The dependence of Rquarki^) on 6 is shown in Fig. 1. The middle curve is the central 
value. The uncertainty of the prediction is described by the shaded region. It has been found 
by adding in squares all the parametric uncertainties mentioned above in points (i)-(iv) and 
in eqns. (0)-(|23). 

3.5 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

8 

Figure 1 . Rquark ('^) R function of 6. 



The ratio Rquarki^) is divergent in the limit which is due to the function fssi^)- In 
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this limit, the physical quantity to consider would be only the sum of b^Xs'y and b^Xg 
decay rates, in which this divergence would cancel out. 

The divergence at 6^1 is very slow. It manifests itself only by a small kink in Fig. 1 
where values of 6 up to 0.999 were included. For 6 = 0.99, the contribution from fssi^) to 
Rquarki.5) IS wcU below 1%. In order to allow easy comparison with previous experimental and 
theoretical publications, we choose 5^°"^ = 0.99 as the first particular value of 6 to consider. 
This value of 6 will be assumed in discussing the "total" B—^Xg'y decay rate below. 

A value of 6 which is closer to what is actually measured is 6 = z = (mc/m;,)^. It 
corresponds to counting only the photons with energies above the charm production threshold 
in the 6-quark rest frame. This is the second value of 6 we will consider below. 

We find 

RquarkiS'^n = (3.12 ± 0.29) X 10-3, (48) 
Rguark{S = z) = (2.86 ± 0.24) X 10"^ (49) 

In both cases, the uncertainty is below 10%. The dominant sources of it are mc/mb and Hb- 
The relative importance of various uncertainties is shown in Table. 1. 





Source 




as{Mz) 








^b,pole 




CKM angles 




2.5% 


1.7% 


6.6% 


5.2% 


0.5% 


1.9% 


2.1% 


6 = z 


2.2% 


1.7% 


3.3% 


6.6% 


0.6% 


1.9% 


2.1% 



Table I. Uncertainties in Rquarki^) due to various sources. 



Let us note that setting 5 = z introduces an additional spurious inaccuracy due to md rrib 
in Rquark- Howcvcr, this additional uncertainty is significantly smaller than the original 
inaccuracy due to mc/mb for 5 uncorrelated with mdmb- 

It is quite surprising that the /ib-dependence is weaker for 5 = z than for 5 = S"^"'^. 
Naively, one would expect larger /Zb-dependence for small values of 6, when the first term 
in eqn. (|32D becomes more and more important. Clearly, S = z is not small enough for this 
effect to show up. The weaker /i;,-dependence ior 6 = z seems to be quite accidental. 



6. In the end, we need to pass from the calculated 6-quark decay rates to the 5-meson decay 
rates. Relying on the Heavy Quark Effective Theory (HQET) calculations we write 

/ xNP xNP\ 

B(5^X,7) = BiB^XceUc) R.^arkiS'^n 1 " ^ + ^ , (50) 



m$ ml 
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where S^i^ and S^J^ parametrize nonperturbative corrections to the semileptonic and radia- 
tive 5-meson decay rates, respectively. 

Following refs. 0, we express 6^^^ and 6^J^ in terms of the HQET parameters Ai and 

5Z = (52) 

where g{z) is the same as in eqn. (^2|) . The poorly known parameter Ai cancels out in the 
r.h.s. of eqn. (|50|) which depends only on the difference 6^1^ — S^J^. The value of A2 is known 
from B*-B mass splitting 

A2 = ^(m|, -m|) ^ 0.12 GeV^ (53) 

The two nonperturbative corrections in eqn. ( |50[ ) are both around 4% in magnitude, 
and they tend to cancel each other. In effect, they sum up to only around 1%. Such a 
small number has to be taken with caution, because the four-quark operators Pi, Pq 
have not been included in the calculation of S^J^. Contributions from these operators could 
potentially give one- or two-percent effects. Nevertheless, it seems reasonable to conclude 
that the total nonperturbative correction to eqn. ( pO] ) is well below 10%, i.e. it is smaller 
than the inaccuracy of the perturbative calculation of Rquarki^^' 



When the above caution is ignored and B{B-^XceVe) = (10.4 ± 0.4)% is used, one 
obtains the following numerical prediction for B-^X^'j branching ratio 

B{B^Xs-f) = (3.28 ± 0.33) x 10"^ (54) 

The central value of this prediction is outside the la experimental error bar in eqn. (0).0 
However, the experimental and theoretical error bars practically touch each other. Therefore, 
we conclude that the present B^Xg'y measurement remains in agreement with the Standard 
Model. This conclusion holds in spite of that the theoretical uncertainty is now more than 
twice smaller than in the previously available leading order prediction [jl|. 

In the measurement of B^X^'y branching fraction, one needs to choose a certain lower 
bound on photon energies. Instead of talking about the "total" decay rate, it is convenient to 



^We identify the 6-quark decay rate given by CLEO with the B-mcson decay rate. Large statistical errors 
in the CLEO result make this identification acceptable. 
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count only the photons with energies above the charm production threshold in the S-meson 
rest frame 

ml -ml ^ 

' 2mB 

Most of the photons in B^Xs'j survive this energy cut, and a huge background from 
charm production is removed. Let us denote the corresponding branching fraction by 

If the 6-quark was infinitely heavy, it would not move inside the 5-meson. The restriction 



(|55|) would then be equivalent to setting 5 = z at the quark level. Since the Fermi motion 



of the 6-quark inside the 5-meson is a 0{A/mb) effect, we can write similarly to eqn. (pO]) 



B{B^X,^r"-' = BiB^X.eu) R,uark{S = z) {l + 0{A/m,)). (56) 

The 0{A/mi,) part of the latter equation can be calculated using specific models for the 
Fermi motion of the 6-quark inside the i?-meson, as it has been done e.g. in ref. [|T^ (see 
also ref. fSll). Performing a similar calculation with use of the NLO values of the Wilson 
coefficients is beyond the scope of the present paper. Such an analysis will be necessary when 
more statistics allows to reduce experimental errors in measurements of B{B—^Xs'^Y^°'"^. 



7. To conclude, we have presented practically complete NLO formulae for b—^Xg^y decay in 
the Standard Model. They include previously published contributions as well as our new 
results for three-loop anomalous dimensions. Our prediction for B—>-Xs'~f branching fraction 
in the Standard Model is (3.28 ± 0.33) x 10~^ which remains in agreement with the CLEO 
measurement at the 2a level. Clearly, an interesting test of the SM will be provided when 
more precise experimental data are available. More importantly, the B—>-Xs'-f mode will have 
more exclusion power for extensions of the SM. 
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Appendix A. 

Here, we give the eight evanescent operators we have used in our anomalous dimension 
computation. Giving them exphcitly is necessary in order to fully specify our renormalization 
scheme, and thus give meaning to the anomalous dimension matrix presented in eqn. (^. 

q 

Pw = (si7^,7M2 7M3 7M47M5 2^"&L)E(^7'^'7'^'7^^7'^^7'^^T'^g)-20P6 + 64P4 

P21 = (sL7Mi7M2 7M3 7M47M5^"cL)(cL7'^i7'^^7'^^7^^7'''T'^&L)-20Pn- 256Pi (57) 

P22 = (si7Mi7M2 7M3 7M47M52^Ci)(cL7'^'7^^7'^^7'^^7^^&L)-20Pi2 - 256P2 

P25 = {sLi,,i,,i,,i,,i,,i^.,iM Y^iqY'Y'Y'Y'Y'Y'Y'q) - 336P5 + 128OP3 

q 

P26 = (sL7Mi7M2 7M3 7M47M57/.e7/.r^"&L) E(97'^^7^^7^^7'^^7^^7'^^7^^T'^g) - 336P, + I28OP4. 



Appendix B. 

Here, we express the functions fij{S) present in eqn. (|32D in terms of the quantities 
Mki{t,u) given in eqn. (27) of ref. Jl^l-Q For each (ij) 7^ (77) and i < j, fij{S) is given by 
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f,,{5) = -(2 - 5ij) E [ du f dt Mkiit, u 



fe,i=i 
k<l 



(58) 



where 



W 



1 
2 
















1 














1 

_ 1 1 

6 2 

16 

-I 8 








1 














1 


2 
1 
1 



used in ref. [17 



P: 



D=A 



k- 



(59) 



In four space-time dimensions, the matrix W transforms our operator basis (H) into the one 



(60) 



k=l 



^ There is a global factor of 2 missing in the expression for Mjs given in eqn. (27) of ref. ||ri 
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